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ABSTRACT 

We compute the cross-correlation between a sample of 14,000 radio-loud AGN 
(RLAGN) with redshifts between 0.4 and 0.8 selected from the Sloan Digital Sky Sur- 
vey and a reference sample of 1.2 million luminous red galaxies in the same redshift 
range. We quantify how the clustering of radio-loud AGN depends on host galaxy mass 
and on radio luminosity. Radio-loud AGN are clustered more strongly on all scales 
than control samples of radio-quiet galaxies with the same stellar masses and redshifts, 
but the differences are largest on scales less than ~ 1 Mpc. In addition, the clustering 
amplitude of the RLAGN varies significantly with radio luminosity on scales less than 
~ 1 Mpc. This proves that the gaseous environment of a galaxy on the scale of its dark 
matter halo, plays a key role in determining not only the probability that a galaxy is 
radio-loud AGN, but also the total luminosity of the radio jet. Next, we compare the 
clustering of radio galaxies with that of radio-loud quasars in the same redshift range. 
Unified models predict that both types of active nuclei should cluster in the same 
way. Our data show that most RLAGN are clustered more strongly than radio-loud 
QSOs, even when the AGN and QSO samples are matched in both black hole mass 
and radio luminosity. Only the most extreme RLAGN and RLQSOs in our sample, 
with radio luminosities in excess of ~ 10 26 W Hz -1 , have similar clustering properties. 
The majority of the strongly evolving RLAGN population at redshifts ~ 0.5 are found 
in different environments to the quasars, and hence must be triggered by a different 
physical mechanism. 

Key words: galaxies: evolution - galaxies: active - radio continuum: galaxies - 
quasars: general 



1 INTRODUCTION 

In recent years, galaxy formation models have become in- 
creasingly interested in the radio AGN phenomenon, be- 
cause it is hypothesized that these objects may regulate the 
star formation history and mass assembly of the most mas- 
sive galaxies and black holes in the Universe. Nearby ra- 
dio galaxies in clusters are observed to inject a significant 
amount of energy into the surrounding gas. As the radio jets 
expand and interact with the surrounding medium, they are 
believed to heat the gas and prevent further accretion onto 
the central galaxy. 

The precise conditions that determine whether an AGN 
develops radio jets/lobes are still a matter of debate. Sev- 
eral studies have shown that the probability for a galaxy to 
become radi o-loud is a strong function of stellar mass and 
redshift (e.g. lBest et~al]|2005Mbonoso et alj|2009h . The role 
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that the environment plays in triggering or regulating the 
RLAGN phenomenon is not as well established. 

iLedlow fc Owenl l| 1996ft found that the fraction of ra- 
dio sources and the shape of the bivariate radio-optical 
was the same for objects in cluster and field environments. 
iBest et al l (|2007t ) found that group and cluster galaxies had 
similar radio properties to field galaxies, but the brightest 
galaxies at the centers of the groups where more likely to 
host radio- loud AGN than other galaxies of the s ame s tellar 
mass. In the local universe, iMandelbaum et al.l (|2009h ana- 
lyzed a large sample of RLAGN at z~0.1. They showed that 
RLAGN inhabit massive dark matter halos (> 10 12 ' 5 M Q ) 
and that, at fixed stellar mass, radio- loud AGN are found 
in more massive dark matter halos than control galaxies 
of the same mass selected without regard to AGN prop- 
erties. This result implies that RLAGN follow a differ- 
ent halo mass - stellar mass relation than normal galax- 
ies. IMandelbaum et al.l (|2009l ) also found that the boost to- 
wards larger halo masses did not depend on radio luminos- 
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itv. lHickox et all (|2009h investigated the clustering in a small 
sample of higher redshift radio-loud AGN selected from the 
AGN and Galaxy Evolution Survey (AGES). They found no 
difference in the clustering amplitude of radio galaxies when 
compared to normal galaxies matched in redshift, luminosity 
and color. 

Most nearby RLAGN lack any of the standard 
accretion-related signatures that would indicate that their 
black holes are growin g significantly at the present day 
l|Hardcastle et al 120061 ). In contrast, quasars are thought 
to be powered by supermassive black holes accreting at 
close to the Eddington rate. Large redshift surveys like the 
Two Degree Field Galaxy Redshift Survey (2dFGRS) and 
the Sloan Digital Sky Survey (SDSS) now provide angu- 
lar positions, accurate photometry and spectra for tens of 
thousands of QSOs. Recent determinations of the quasar 
two-point correlation function have demonst rated that at 
z < 2 . 5 quasars cluster like normal L* galaxies l|Croom et al.l 
I2OO5I ; ICoil et all [2007h and populate dark matter halos of 
~ 10 12 Mq , with the clustering only weakl y dependent on lu - 
minosity, color and virial black hole mass (|Shen et alJl2009h . 

As one moves out in redshift, the radio-loud AGN popu- 
lation evolves very rapidly in radio luminosity. Whether the 
RLAGN population also evolves strongly in black hole accre- 
tion rate, is considerably less clear. In particular, our under- 
standing of whether there is a relationship between powerful, 
high redshift radio-loud AGN and quasars is quite sketchy. 
Around 10% of the quasar population is radio-loud. Numer- 
ous investigations have found that radio-loud quasars and at 
least some powerful radio galaxies share a number of com- 
mon characteristics, such as excess infrared emission, com- 
parable radio morphologies and luminosities, optical emis- 
sion lines, large evolutionary rates, and host galaxies with 
similar properties. It has thus been tempting to link both 
phenomena under the hypothesis that they are t he same 
active nuclei viewed at diffe rent orientations (e.g. iBarthell 
1 19891 ; lUrrv fc Padovanilll995l ) . 

A few facts are believed to be key in any attempt 
to understand the transition from the population of low- 
luminosity radio AGN produced by weakly accreting black 
holes at low redshifts, to a population of high-luminosity 
radio AGN that may be p roduced by strongly ac creting 
black holes at high redshifts. iFanaroff fc Rilevl (|l974T ) found 
an important correlation between radio morphology and 
radio power: low luminosity sources (Fanaroff-Riley Class 
I, FRI) show emission peaking close to the nuclei that 
fades toward the edges, whereas more luminous sources 
(Fanar off-Riley Class I I, FR II) are brightest toward the 
edges. iHine fc Longairl (|l979h discovered that radio galax- 
ies could also be classified according to the strength of 
their optical emission lines: low excitation (weak-lined) ra- 
dio galaxies or LERGs, and high excitation (strong-lined) 
objects or HERGs. Modern unification models usually as- 
sociate quasars with the most powerful HERGs, and low 
luminosity LERGs with BL Lac objects. Although there is 
a notable correspondence between RLAGN luminosity, mor- 
phology and spectral type, i.e. lower luminosity FRIs with 
LERGs, and higher luminosity FRIIs with HERGs, the cor- 
relations between these properties are not straightforward. 
There are populations of FRI sources with high excitation 
nuclear lines, and conversely, FRII galaxies with low excita- 
tion spectra are also common. 



It has been known for years that very high redshift (z > 
2), powerful radio galaxies are often surro unded by galaxy 
overdensities with sizes of a few Mpc (e.g. iPentericci et al.l 
120001 : iMilev et al.ll2006r ). Since we know that quasars at the 
same redshift are clustered like normal L, galaxies, this 
would seem to throw some doubt on a simple unified scheme 
for explaining both phenomena. 

In view of this highly complex situation, a more statis- 
tical approach to comparing the properties of quasars and 
radio galaxies may yield further insight. In this paper we 
present measurements of the projected cross-correlation be- 
tween a sample of 14,000 radio-loud AGN with a median 
redshift of z = 0.55 with the surrounding population of 
massive galaxies (M, > 10 11 Mq). The large size of our 
samples allows us to investigate in detail how clustering de- 
pends on stellar mass and on radio luminosity. By compar- 
ing the RLAGN clustering with results from control samples 
matched in redshift, luminosity and mass, we isolate the ef- 
fect that the radio AGN phenomenon has on the cluster- 
ing signal. We cross-correlate radio quasars drawn from the 
SDSS with the same reference sample of massive galaxies. 
Again, by using control samples matched in black hole mass 
and radio luminosity, we ensure that we compare RLAGN 
and RLQSOs in as uniform a way as possible. 

This paper is organized as follows. In Section 2 we de- 
scribe the surveys and samples used in this work. In Section 
3 we explain the methodology adopted to calculate the two- 
point correlation function. Section 4 presents the results on 
radio-loud AGN and quasar clustering. Finally, in Section 
5 we summarize our results and discuss the implications of 
this work. 

Throughout the paper we assume a flat ACDM cosmol- 
ogy, with Q m ~ 0.3 and SIa = 0.7. Unless otherwise stated, 
we adopt h = Ho /(100 km s _1 ) and present the results in 
units of Mpc with h = 1. 

2 DATA 

2.1 The MegaZ-LRG Galaxy Catalogue 

The Sloan D i gital Sky Survey (|York et all |2000| : 
IStoughton et al.l 120021 ') is a five-band photometric and 
spectroscopic survey that has mapped almost a quarter of 
the whole sky, providing precise photometry for more than 
200 million objects and accurate redshifts f or about a mil- 
lion g alaxies and quasars. The MegaZ-LRG |Collister et al.l 
20Q7S) is a photometric redshift catalogue based on imaging 
data from the fourth Data Release (DR4) of the SDSS. 
It consists of ~1.2 million Luminous Red Galaxies (LRG) 
with limiting magnitude i < 20 over the redshift range 
0.4 < 2 < 0.8. MegaZ adopts various color and magnitude 
cuts to isolate red galaxies at 0.4 < z < 0.8. The cuts 
are very similar to those adopted by the '2dF-SDS S 
LRG and Quasar' project (2SLAQ. ICannon et all l2006h . 
Accurate photometric redshifts are available for the entire 
LRG sample. These are derived using a neural network 
photo metric redshift estimator (ANNz, Collis ter fc Lahavl 
120041 ) that was trained using a sample of ~ 13000 LRG with 
spectroscopic redshifts selected from 2SLAQ. The r.m.s. 
average photometric redshift error for all the galaxies in the 
sample is a rms = 0.049. 
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2.2 The Radio- loud Galaxy Sample 

By combining the optical MegaZ-LRG catalo gue with data 
from the NRAO VLA Sky Survey (NVSS; Icondon et aD 
1998) and the VLA Fain t Images of the Radi o Sky at Twenty 
Centim eters (FIRST; Beck er et al.l Il995l ). iDonoso et al.l 
(2009) constructed a catalogue of 14453 radio-loud AGN 
with 1.4 GHz fluxes above 3.5 mjy. The cross-matching 
method utilized a collapsing algorithm to identify multiple- 
component FIRST and NVSS sources and the method was 
optimized to take advantage of both surveys. NVSS has suf- 
ficient surface brightness sensitivity to provide accurate flux 
measurements of extended radio sources with lobes and jets. 
On the other hand, the superior angular resolution of FIRST 
is crucial to identify the central core component of each ra- 
dio source and to provide a robust association between the 
radio source and the optically-identified host galaxy. 

Monte-Carlo simulations were used to estimate the reli- 
ability (~98.3%) and completeness level (95%) of the cata- 
logue. The vast majority of the detected radio AGN (78.6%) 
are single component sources in both NVSS and FIRST. 
There is, however, a significant fraction of objects without 
(catalogued) high S/N FIRST detections (~8%), so the au- 
thors introduced a method for analyzing radio maps that 
allowed them to dig deeper into the FIRST survey and to 
use lower S/N detections to pinpoint the location of the 
host galaxy. We ref er the reader to the original paper by 
IDonoso et al.l (|2009) for a detailed description of these pro- 
cedures and the matching algorithm. 



2.3 The Radio- loud and Radio-quiet Quasar 
Samples 

In this work we use quasars selected from the fourth edi- 
tion of the spectroscopic quasar catalogue ijSchneider et al.l 
I2007D . This contains 77,429 quasars drawn from SDSS DR5, 
with luminosities larger than Mi = —22, that have at least 
one broad emission line with FWHMM000 km s~ x in their 
spectra. The catalogue also identified radio-loud quasars 
with FIRST components within a 2 arcsec radius. 

Most of the objects target ed as quasars w e re init ially se- 
lected using the algorithm of I Richards et all l|2002l ). which 
pick candidates using ugriz broadband photometry and by 
matching with unresolved FIRST sources. As the survey pro- 
gressed, the quasar selection software was modified to im- 
prove its efficiency at high redshift. This is reflected in two 
spectroscopic target selection flags listed as TARGET and 
BEST (for the final algorithm). Photometry of quasars is 
also available in two versions, TARGET measurements (val- 
ues used at the time of targeting) , and BEST measurements 
(values derived with the latest pipeline). We note that the 
selection of UV-excess quasars at low redshifts [z < 3) has 
remained essentially unchanged, so that only small differ- 
ences arise from using TARGET or BEST versions. In addi- 
tion, the bias introduced by selection of targets via FIRST 
radio detections is significant only at high redshift. 

In this work we are interested in cross-correlating 
the quasars with the LRGs described above. We there- 
fore selected an homogeneous quasar sample consisting of 
all quasars with 0.35 < z < 0.78 and psf magnitudes 
in the range 15 < i < 19.1. We only consider primary 
objects (primary = 1) with point source morphology 
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Figure 1. Top: normalized distribution of radio luminosity 
(Pl.4GHz) corresponding to radio-loud AGNs (red), and to QSOs 
detected down to 1 mjy in the FIRST survey (blue). The vertical 
line at 10 25 W Hz^ 1 marks the adopted threshold between radio- 
quiet and radio-loud QSOs (dotted). Also shown is the ratio of 
the number of radio-loud quasars relative to radio AGN (scale on 
the right axis). Bottom: distribution of black hole mass (Mj^) for 
radio-loud AGNs and radio-loud QSOs. 



(morphology = 0), that were also targeted as primary sci- 
ence objects (scienceprimary = 1). This yields a sample of 
7128 quasars. 

Of these 7128 quasars, 684 (9.6%) have radio identifica- 
tions in the FIRST survey down to the 1 mjy flux density 
limit. One issue that could affect the derived radio lumi- 
nosities of the QSOs in our sample is that a fraction of 
them present a truly extended FRII-like morphology, and 
the total radio flux is distributed over many components. 
The exclusion of such structures might lead to an underesti- 
mation of the total radio luminosity. We visually examined 
NVSS/FIRST radio maps of the 678 QSO with FIRST de- 
tections and added the NVSS fluxes of the associated com- 
ponent^), if present, or of the FIRST component(s) when 
no NVSS source was found in the nearby. For some radio 
QSO the derived radio luminosities increase by a factor of 
~2-3. Nevertheless, we note that we repeated the clustering 
analysis described in Section ^. 2l using only the central (core) 
component flux, and we verified that this has no significant 
influence on any of our results. 

According to convention, we define radio-loud quasars 
as those with total integrated 1.4 GHz radio power (af- 
ter adding all associated components) above 10 25 W Hz -1 . 
With this definition, there are 307 radio-loud quasars in our 
sample. We consider objects below this luminosity (or non- 
detections) as radio-quiet quasars. 
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2.4 Sample Properties 

Sh en et al.l (|200Sl } have derived virial black hole mass es- 
timates for SDSS DR5 quasars. These are based on H/3, 
Mgll, and CIV emission lines, and the continuum lumi- 
nosities around these lines. We adopt these estimates for 
our quasar sample (at z < 0.7, these are mostly derived 
from H/3) . For RLAGN we adopted the relation betwee n 
black hole mass-bulge mass derived bv lHaring fc Rbd (|2004 ). 

= 0.0014Mf, u ; 9e , where we replace Mf, u i gf , by the stellar 
mass of the galaxy. At the lower end of our galaxy mass dis- 
tribution (~ W 11 Mq), use of the stellar mass instead of the 
bulge mass may cause the black hole mass to be overesti- 
mated by a factor of ~ 1.2 — 1.4. We note that the majority 
of RLAGN in our sample are more massive than this. 

For reference, Figure [1] shows the radio luminosity and 
black hole mass distributions derived for all the radio-loud 
AGN and radio-loud QSO in our samples. In the upper 
panel, we also plot the ratio of the number of RLQSOs 
to RLAGN and show that this increases from ~ 1% at 
10 25 W Hz" 1 up to ~ 50% at 1 27 W Hz" 1 . Th is is 
broadly consistent with the results of lLawrenc 1 (|l99lh who 
found similar relative proportions between broad-lined and 
narrow-lined 3CR sources. One scenario that has been intro- 
duced to explain this varying fraction (at least at luminosi- 
ties where the relative numbers of RLQSO and RLAGN ar e 
similar) is the receeding torus model (e.g. ISimpsonlll998l ). 
in which the inner radius of the obscuring torus (which is 
identified with the dust sublimation radius) scales with lu- 
minosity as L 5 . This model predicts that a larger fraction 
of more luminous objects are classified as quasars. 



3 CLUSTERING ANALYSIS 

3.1 The Cross-correlation Function 

A standard way to characterize the clustering of galaxies is 
with the two-point correlation function which measures 
the excess in the numbers of pairs with separation r in a 
volume dV, respect to a random di stribution wit h the same 
mean number density of objects n (Peebles 1980). This can 
be expressed as 



dP = n[l + £(r)]dV 2 



(1) 



Objects are said to be clustered if £ > 0. The amplitude and 
shape of the correlation function yield a variety of different 
information. On scales larger than a few Mpc, the ampli- 
tude is a measure of the mass of dark matte r halos that the 
galaxies inhabit (e.g. lSheth fc Tormen| [l999). On intermedi- 
ate scales, the shape of the correlation function is sensitive 
to how galaxies are distributed within their halos (|Li et al.l 
l2006bT ). while at scales smaller than a few hundre d kpc it 
prob es processes such as mergers or interactions (|Li et al.l 
l200Sh . 

Several estimators for the (auto)correlation function 
have been proposed in the literature. In this work we cal- 
culate the a uto-correlation fu nction of the LRGs using the 
estimator of iHamiltonl J 19931 ). 



where DD(r), RR(r), DR(r) refers to the normalized num- 
ber of (LRG-LRG), (random-random), and (LRG-random) 
pairs as function of the spatial separation r (see the next sec- 
tion for details about the construction of the random sam- 
pie). 

To estimate the cross-correlation function of radio-loud 
AGN or quasars with the MegaZ-LRG galaxy sample, we 
count the number of LRGs around each AGN or quasar as 
a function of distance, and divide by the expected number 
of pairs for a random distribution, 



CD(r) 
CR(r) 



1 



(3) 



where CD(r) stands for the number of (RLAGN/QSO- 
LRG) pairs, CR(r) is the number of (RLAGN/QSO- 
random) pairs, and the quantities have been normalized by 
the number of objects in the LRG and random catalogues. 
The advantage of our procedure is that it does not require 
full knowledge of the QSO or RLAGN selection function. 
Only the LRG selection function is needed for the construc- 
tion of the random sample, and this is well quantified. An- 
other reason for calculating cross-correlations rather than 
auto-correlations, is that it allows us to overcome shot noise 
when the sample size is small. We note that the LRG sam- 
ple (D in the notation above) remains fixed throughout this 
work. The error bars of the auto and cross-correlation func- 
tions are calculated via statistical bootstrapping by drawing 
n — 100 random samples with replacement. 

In practice, photometric redshift errors as well as dis- 
tortions due to peculiar velocities along the line of sight will 
introduce systematic effects in our estimate of {(r). There- 
fore, to recover real-space clustering properties we decom- 
pose £ in two directions, along the line of sight (71") and 
perpendicular to it (r p ). Integrating over the 7r-direction al- 
lows to define the projected two-point cross-correlation func- 
tion w v (r v ), a quantity that is independent of such distor- 
tions (|Davis fc Peebleslll983h. A detailed description of the 



method can be found in iLi et al.l (|2006bh 



DD(r)RR(r) _ 
W ' [DR(r)] 2 



(2) 



3.2 Construction of the Random Sample 

The random sample used in estimating the cross-correlation 
function should have the same selection effects as the ob - 
served galaxies. We follow the method bv lLi et al.l (|2006al ) : 
we take observed LRG sample inside the coverage mask of 
SDSS DR4 and randomly re-assign the sky coordinates of 
each galaxy. All other quantities such as like redshift, stel- 
lar mass and luminosity are kept fixed. Because the survey 
covers a very wide area (>6000 deg 2 for SDSS DR4), this 
procedure is sufficient to remove any coherence in the radial 
direction and it ensures that the geometry of the random 
catalogues are exactly the same as the real one, and that all 
redshift-dependent selection effects are accounted for. We 
generate N=10 random samples in this way. 



4 RESULTS 

4.1 Radio-loud AGN clustering 

It is well known that the clustering amplitude of galaxies 
varies as a function of mass, luminosity and redshift. Radio 
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Figure 2. Normalized distributions of i-band absolute magni- 
tude, redshift and stellar mass for radio-quiet LRGs (dotted), 
radio-loud AGN (histogram), and control radio-quiet LRGs (large 
dots). 

AGN are usually hosted by very ma ssive, > 3L* galaxies 
l|Best et al.ll2005l ; iDonoso e't all 120091 ). To take this into ac- 
count, we select control samples of radio-quiet MegaZ-LRG 
galaxies with redshifts, stellar masses and absolute magni- 
tudes that closely match the radio AGN sample. For each 
RLAGN we randomly select 10 radio-quiet LRG (or 5, de- 
pending on the number of available candidates) within a 
tolerance of Az = 0.02 in redshift, AM = 0.1 in log stel- 
lar mass, and AM, = 0.05 in absolute magnitude, where 
Mi is the extinction and k-corrected i-band absolute magni- 
tude. Figure [2] shows the distributions of these parameters 
for radio-loud, radio-quiet and control objects. 

Using the methods described in the previous section, we 
first calculate the auto-correlation function for our reference 
sample of luminous red galaxies. We then cross-correlate the 
radio-loud AGN with the LRG parent sample. This is shown 
in Figure [3] where it can easily be appreciated that RLAGN 
are significantly more clustered than the LRG population 
on all spatial scales. The two terms of the clustering signal, 
corresponding to galaxies within the same halo (1-halo) and 
in different halos (2-halo), are clearly visible with the tran- 
sition occurring around 1 Mpc The boost in clustering 
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Figure 3. Projected cross-correlation function w p (r p ) between 
radio-loud AGN and MegaZ luminous red galaxies (red, solid) in 
the range 0.3 to 30 Mpc h~ 1 . The LRG-LRG auto-correlation 
is indicated by small diamond symbols. Also shown is the cross- 
correlation of a control sample of radio-quiet LRG (red, dashed) 
with the same distribution of redshifts, luminosities and stellar 
masses as the radio-loud population. The bottom panel shows the 
ratio of w p (r p ) for the RLAGN to that for the control sample. 

signal exhibited by the RLAGN is significantly stronger on 
scales less than 1 Mpc hT 1 , which tells us that RLAGN must 
occupy special positions within their dark matter halos. We 
intend to model this in more detail in upcoming work. 

If we compare the clustering of RLAGN with that of 
control galaxies with the same redshifts, luminosities and 
stellar masses, we see that RLAGN are still significantly 
more clustered. The ratio between the cross-correlations 
w p (r p ) corresponding to RLAGN and its corresponding con- 
trol radio-quiet sample, is plotted in the bottom panel. This 
proves that the probability of a galaxy to become radio-loud 
depends on environment as well as on black hole or galaxy 
mass. 

4-1.1 Dependence on Stellar Mass 

We split the RLAGN sample into two subsamples with 
log(M/M Q ) < 11.3 and with log(M/M Q ) > 11.6. We also 
applied the same split to the corresponding control sam- 
ples. The resulting cross-correlations are plotted in Fig- 
ure 2] As expected, more massive radio galaxies are more 
strongly clustered on all scales. When compared to control 
galaxies, both subsamples show roughly the same relative 
clustering strength on scales larger than 1-2 Mpc h . On 
small scales the difference between control and data sam- 
ples is more significantly boosted for RLAGN in less massive 



6 Donoso et al. 




control all RLAGN 
control log(M/M aun ) < 11.3 
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r p / Mpc h" 1 

Figure 4. Top: projected cross-correlation function w p (r p ) be- 
tween radio-loud AGN and MegaZ luminous red galaxies (red) in 
the range 0.3 to 30 Mpc h- 1 . The green and blue lines indicate 
the cross-correlation of massive objects with log(M/M0) > 11.6, 
and of less massive systems with log(M/M©) < 11.3. The LRG- 
LRG auto-correlation function is shown for reference (diamond 
symbols). Middle: cross-correlation of control samples of radio- 
quiet LRGs that have the same distribution of redshift and stellar 
mass as the radio-loud systems. Bottom: ratio of w p (r p ) between 
RLAGN and their corresponding control samples. Note the curves 
are slightly shifted along the x-axis to improve the visibility. 

galaxies. These results ar e in good agreement with those of 
iMandelbaum et al l (|2009T l for RLAGN at lower redshifts. 

We now investigate how the clustering of RLAGN and 
their control galaxies varies as a function of stellar mass. 
We fit two power laws of the form w(r p ) = A r^ 1 7 ' to 
the cross-correlation function, one over the range 0.1 < 
r p < 0.8 Mpc and the other over the range 1 < r p < 
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Figure 5. The dependence of the cross-correlation amplitude of 
RLAGN and control galaxies on stellar mass. The amplitude is 
computed by fitting a power law with fixed exponent (see text for 
details). Results are shown for RLAGN (solid lines) and their 
corresponding control radio-quiet LRG (dotted, dashed lines). 
Fits are calculated at two different spatial scales, 0.1 < r p < 
0.8 Mpc h' 1 (blue) and 1 < r p < 20 Mpc h' 1 (red). 

20 Mpc ft -1 . This division allows us to quantify separately 
the clustering signal contributed by LRGs within the same 
halo as the RLAGN and by LRGs residing in different halos. 
For the complete RLAGN sample, the best fitting parame- 
ters are A = 233.9 ± 15 and 7 = 2.18 ± 0.05 on scales less 
than 1 Mpc h' 1 and A = 173.2 ± 10 and 7 = 1.81 ± 0.05 on 
larger scales. We then divide the sample into 8 mass bins and 
perform new fits, keeping the slope of the power law fixed 
and allowing the normalization to vary. Figure [5] shows the 
cross-correlation amplitudes as a function of stellar mass for 
RLAGN and the radio-quiet control sample. As can be seen, 
the ratio between the clustering amplitude of the RLAGN 
and the control galaxies depends both on stellar mass and on 
the scale at which the clustering is evaluated. On scales less 
than 1 Mpc h^ 1 , there is a relatively strong dependence of 
the ratio on stellar mass, with RLAGN in low mass galaxies 
clustered much more strongly than the controls, but RLAGN 
in high mass galaxies clustered similarly to the controls. On 
larger scales, there is a much weaker trend with mass. 

4-1-2 Dependence on Radio Luminosity 

We now investigate if there is any dependence of 
RLAGN clustering o n the luminosity of the radio source. 
iPrestage fc Peacock! (119881 ) studied the local galaxy density 
around radio galaxies at z < 0.25, finding that weak FRI 
sources are typically found in de nser regions c ompa red to 
the more luminous FRII sources. I Yates et al.l (|l989T ) (and 
later iHill fc Lillvlll99ll) extended such studies to higher red- 
shifts, concluding that powerful radio galaxies at z ~ 0.5 
are typically found in environments three times richer than 
their counterparts at z ~ 0.2, but also that the most lu- 
minous objects (Pi78MHz ~ 10 27 ' 1 W Hz -1 sr _1 ) occupy 
richer environments than the weaker objects (Pi78MHz ~ 
10 26 1 W Hz" 1 sr" 1 ). However, given the limitations of the 
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samples available, they were unable to determine if such 
clustering trends were primarily dependent on r edshif t or on 
radio luminosity, or on a combination of both. iBestl (|2004l ) 
studied the density of galaxies around nearby radio-loud 
AGN. He found a positive correlation between local density 
and radio luminosity for RLAGN without emission lines, 
but found that RLAGN with emission-lines tended to avoid 
regions of high density. 

In this study, we split our RLAGN sample into a low- 
luminosity subsample with log(Pi.4GHz [W Hz - ])< 24.6, 
and a high-luminosity subsample log(Pi.4GHz [W Hz _1 ])> 
25.7. These cuts allow us to sample the faint and bright end 
of the radio luminosity distribution. We again build control 
samples in the same way as before and we present the cross- 
correlation results in Figure [U The top panel of Figure [5] 
shows that low luminosity RLAGN are more clustered than 
high luminosity systems at all scales. When compared to 
control samples, this "boost" in clustering is only visible at 
scales larger than ~1 Mpc /i -1 . 

To quantify the variation of clustering with radio lu- 
minosity in more detail, we once again proceed by fitting 
a power law to the cross-correlation functions for RLAGN 
subsamples split by radio luminosity. We fit separate power 
laws on scales below and above 1 Mpc h~ . The varia- 
tion in the clustering amplitude with luminosity is plotted 
Figure [7] Two interesting features can be observed. First, 
the clustering amplitude of radio galaxies on large scales is 
only very weakly anti-correlated with radio power. On small 
scales, the clustering increases with radio luminosity, peaks 
at log(Pi.4GHz[W Hz~ 1 ])~ 25.3, and then decreases for most 
lum inous radio sources. 

iBarthel fc Arnaudl l|l99fj) argue that the confining ef- 
fect of a dense intracluster medium reduces the adiabatic 
loses of radio lobes, leading to higher levels of synchrotron 
emission. Thus, a dense environment may provide a more 
effective 'working surface' for the lobes, giving rise to the 
positive correlation between small-scale clustering ampli- 
tude and radio luminosity observed in Figure [7] for sources 
with log(Pi.4GHz[W Hz~ 1 ])< 25.3. Alternatively, higher ra- 
dio luminosities in denser environments may be a result 
of increased jet powers resulting from the higher cooling 
rates in these denser regions. Why does the clustering am- 
plitude drop for radio sourceswith luminosities higher than 
this value? As we will argue in the next section, a radio 
luminosity of log(Pi.4GHz [W Hz _1 ])~ 25.3 may mark the 
beginning of a transition to a population of AGN that are 
more similar to the quasars, which as we will show, are sig- 
nificantly less clustered that the RLAGN. 

4.2 Quasar clustering and AGN Unification 

In this section, we compare the clustering of radio galaxies 
and quasars at z~0.5. Our goal is to develop a better un- 
derstanding of the relationship between these two types of 
active galaxy. 

AGN unification models provide an appealing way to 
account for the diversity of the observed AGN population. 
The basic hypothesis is that the observed characteristics of 
AGN depend mainly on their orientation relative to the line- 
of-sight. Comprehensive reviews of unifi cat ion models can be 
found in lBarthel|[l989l . lAntonuccilll993l or lUrrv fc Padovanil 
1 19951 . 




control all RLAGN 
control log(P) < 24.6 
control log(P) > 25.7 
LRG auto 




r p / Mpc IV 1 

Figure 6. Top: projected cross-correlation function w p (r p ) 
between radio-loud AGN and MegaZ luminous red galaxies 
(red) in the range 0.3 to 30 Mpc h~ x . Green and blue 
lines indicate the cross-correlation of luminous objects with 
log(Pi.4GHz[W Hz~ 1 ])> 25.7, and of less powerful AGN with 
log(Pi. 4G Hz[W Hz~ 1 ])< 24.6. The LRG-LRG auto-correlation is 
shown for reference (diamond symbols). Middle: cross-correlation 
of control samples of radio-quiet LRGs that have the same dis- 
tribution of redshift and stellar mass as the radio-loud systems. 
Bottom: ratio of w p (r p ) between RLAGN and their correspond- 
ing control samples. Note the curves are slightly shifted along the 
x-axis to improve the visibility. 

In this paper, we attempt to test one fundamental re- 
quirement of the unification scheme of radio-loud objects, 
namely that the environment of radio galaxies and radio 
quasars should be statistically identical. We note that pre- 
vious work has already suggested that low excitation ra- 
dio galaxies (which include most FRI sources, but also 
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Figure 7. The dependence of the cross-correlation amplitude on 
radio luminosity. Results are shown for both, RLAGN (solid lines) 
and their corresponding control radio-quiet LRGs (dotted, dashed 
lines). Fits are calculated for two different ranges in scale: 0.1 < 
r p < 0.8 Mpc h- 1 (blue) and 1 < r p < 20 Mpc h' 1 (red). 



a significant fraction low luminosity FRII radio galaxies) 
do not participate in the same unificati on framework as 
quasars or broad line radio galaxies (e.g. iHardcastlel 1 20041; 
lHardcastle et al1l2007t ). We will therefore confine our atten- 
tion to the most luminous radio-loud galaxies and radio-loud 
quasars in our sample, i.e. both with luminosities in excess 
of 10 25 WHz" 1 . 

Up to now, observational evidence has not yielded con- 
clusive evidence as to whether powerful RLAGN and RLQ- 
SOs cluster in the same way. The first problem is that the 
available samples have been small. In the local universe, 
powerful radio galaxies with log(Pi.4GHz [W Hz- 1 ])- 26 have 
typical comoving densities of 10 -8 Mpc -3 dex -1 at z~0.1, 
so large vol umes are required t o detec t a significant number 
of sources. ISmith fc Heckmanl l|l990fl studied the environ- 
ments ~30 low redshift radio quasars and powerful radio 
galaxies, concluding that both populations were clustered 
in much the same way as radio-quie t QSO s. At higher red- 
shifts (0.3 < z < 0.5). lYates et all (|l989l ) also found that 
the environments of radio galaxies and radio-loud quasars 
were simi lar, with higher l uminosity systems slightly more 
clustered. iBarr et all l)2003h found that luminous radio-loud 
quasars exist in a variety of environments including rich clus- 
ters, compact groups and in low-density environments. 

In this work we calculate the cross-correlation func- 
tion between radio-quiet and radio-loud quasars, and the 
same reference sample of LRGs used in Section 4.1. The 
resulting w p (r p ) are plotted in Figure [8] As can be seen, 
there is no significant difference in clustering strength be- 
tween radio-loud and radio-quiet quasars. It is interesting 
that the clustering strength of RLAGN seems to be larger 
than that of radio-loud quasars on all scales, and particu- 
larly at r p <1 Mpc ft -1 . The mean relative bias of radio- loud 
quasars respect to the LRG population remains roughly con- 
stant at ~0.7, while the bias of RLAGN varies strongly from 
~2.5 to ~1.5 below 1 Mpc h~ , and then stays relatively 
constant at larger scales. 




0.1 1.0 10.0 

r p / Mpc h 1 



Figure 8. Projected cross-correlation function w p (r p ) between 
quasars and LRGs (dotted, black), and between radio- loud 
quasars and LRGs (blue). For comparison, we plot the cross- 
correlation of radio-loud AGN and LRGs (red), as well as 
the LRG-LRG auto-correlation (small diamonds). The grey 
sh aded area indicates t he Q SO-LRG cross-correlation derived 
bv lPad manabha n et al.l ll2008r i. The bottom panel shows the ra- 
tio of iv p (r p ) respect to the LRG auto-correlation. The anal- 
ysis is restricted to sources with integrated luminosities above 
10 25 W Hz" 1 . 



We note that our quasar/LRG cross-correlation 
function agrees ex t remel y well with that derived by 
iPadmanabhan et al.1 (|2008l ). IShen et all l|2009l) analyzed the 
clustering of radio-loud and radio-quiet quasars in SDSS 
DR5 at 0.4 < z < 2.5 and found that radio quasars cluster 
more strongly than radio-quiet quasars with the same black 
hole masses. As we will show in section 4.2.1., matching the 
radio-quiet and radio-loud quasar sample in black hole mass 
does not alter our conclusio n. We speculate that disagree- 
ment with Sh en et al] (|2009l ) may arise because we consider 
a muc h narrower range in redshift. We note that Wold et al.l 
(2000) also found little difference between the environments 
of radio-loud and radio-quiet quasars over roughly the same 
redshift range as that probed in this study. 

We conclude, therefore, that powerful radio galaxies ap- 
pear to be hosted by very massive halos, more massive than 
their quasar counterparts. In principle, this suggests that 
the unification scheme for the two classes of AGN is not as 
straightforward as first thought, and additional parameters 
other than orientation are required to explain the difference 
between RLAGN and RLQSOs. 
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1.0 

r p / Mpc h" 1 

Figure 9. Projected cross-correlation function w p (r p ) between 
radio-loud quasars and LRGs (blue). Also shown is the cross- 
correlation of a control sample of radio-loud AGN (green) selected 
to have a similar distribution of M^/Pi^ghz as the radio-loud 
quasars. For comparison we plot the cross-correlation of RLAGN 
and LRG (red), and the LRG-LRG auto-correlation (small dia- 
monds). The bottom panel shows the ratio of w p (r p ) to the LRG 
auto-correlation. The analysis is restricted to sources with inte- 
grated luminosities above 10 25 W Hz -1 . 



4.2.1 Black Hole Mass 

One such parameter could be the mass of the black hole. 
Some observational evidence supports the idea that radio 
jet power might be closely related to the mass of the black 
hole and its accretion rate. Links between radio luminos- 
ity and black hole mass hav e been found in r adio galax- 
ies (iFranceschini et alJ [l99§) , and in quasars (|Lacv et al.l 
l200ll ; lBorosonll2002| y However, ot her authors have argued 
against such strong correlations (iHol |2002| ; ISnellen et al.l 



120031 ; iMetcalf fc Magliocchettill2006 

To control for the effect of black hole mass, we con- 
structed a sample of RLAGN with a similar distribution in 
Mbh/Pi.4GHz as that of radio-loud quasars. The parameter 
Mb/j/Pi.4GHz can be considered as a kind of inverse Edding- 
ton ratio that measures how much radio emission per unit 
black hole mass is produced by the jet. Figure [9] shows the 
resulting cross-correlations. A slight decrease in w p (r p ) is 
observable at scales above 1 Mpc ft -1 for the RLQSO sam- 
ple, but the effect is of low significance. Figure [lOl shows the 
distribution in Mf, h /Pi.4GHz of the RLAGN and RLQSO 
samples before and after the matching procedure. 

We find that the difference in clustering does not change 
for the samples that are matched in Mbh/PlAGHz, meaning 




Figure 10. Normalized distributions of M;,h/Pi.4GHz for radio- 
loud AGN (dotted), radio-loud QSOs (histogram), and control 
radio AGN (large dots) selected to have a similar distribution in 
Mt,h/Pi.4GHz as the radio-loud quasars. 



that black hole mass has a negligible influence in driving the 
observed differences between the clustering of radio quasars 
and radio galaxies. From Figure [1] it can be seen that some 
radio-loud quasars are hosted by black holes more massive 
thanlO 9 ' 3 M , which are not present in the RLAGN popula- 
tion (we suspect that errors in the virial black hole mass esti- 
mates are to blame) . We have repeated the cross-correlation 
analysis of radio quasars with black hole masses in the range 
10 8 M Q < M bh < 10 9 M and find that this makes no dif- 
ference to our results. 



4-2.2 Radio Luminosity 

It is also interesting to investigate whether clustering dif- 
ferences between radio-loud AGN and quasars depend on 
radio luminosity. To test this, we build control samples us- 
ing the same methodology as before, but this time matching 
in log(Pi.4GHzj3- Figures [TT1 and [T2l show the correspond- 
ing cross-correlation functions and radio luminosity distri- 
butions of the matched samples. The clustering of RLAGN 
remains essentially unchanged. 

We now calculate cross-correlation functions for 
radio-loud quasars of increasing radio luminosities 
[log(Pi.4GHz)>25.0,>25.5,>25.75]. We do find an in- 
crease in clustering strength as function of radio power on 
all scales in the range 0.1 < r p < 20 Mpc h^ 1 . This is plot- 
ted in Figure 1131 where we compare the cross-correlation 
amplitude of radio-loud AGN and QSO. The amplitude is 
calculated using a single power-law fit over the entire range, 
since the correlation function of quasars does not exhibit a 
clear break at a scale of ~1 Mpc h~ , as is the case for radio 
galaxies. We find that RLAGN are more strongly clustered 
than RLQSO at all radio luminosities that we are able 



1 We note that a fraction of the RLQSOs will be core-dominated, 
so that a fraction of the luminosity of some sources will be due 
to beaming. This would affect the matching in radio luminosity 
between beamed and non-beamed objects. However, because of 
the weak dependence of clustering amplitude on radio luminosity 
for the quasars, this effect does not influence our conclusions 
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r p / Mpc h" 1 

Figure 11. Projected cross-correlation function w p (r p ) between 
radio quasars and LRG (blue). Also shown is the cross-correlation 
of a control sample of radio-loud LRG (green) selected to have 
a similar distribution of log(Pi.4GHz) as in radio quasars. For 
comparison we plot again cross-correlation of RLAGN and LRG 
(red), and the LRG-LRG auto-correlation (small diamonds). The 
bottom panel shows the ratio of w p (r p ) respect to the LRG 
auto-correlation. The analysis is restricted to sources with in- 
tegrated luminosities (after adding all associated components) 
above 10 25 W Hz" 1 . 
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Figure 12. Normalized distributions of log(P! .4GHz) f° r radio- 
loud AGN (dotted), radio-loud QSO (histogram), and control ra- 
dio AGN (large dots) selected to have a similar distribution as of 
radio-loud quasars. 



350 - 



300 ■ 



200 ■ 



150 



50 ■ 



t 1 1 1 1 r- 

O O RLAGN 0.1 

□ ■ ■ ■□ RLAGN 1 
<S O RLQSO 0.1 

□ RLQSO 1 
RLQSO 0.1 



< 0.8 Mpc 

< 20 Mpc 

, < 0.8 Mpc 
, < 20 Mpc 
, < 20 Mpc 




24.5 



25.0 25.5 
39,o(L,.4ghz / W Hz"') 



26.0 



Figure 13. Change of the cross-correlation amplitude for 
RLAGN (red) and RLQSO (blue), obtained by fitting a 
power law with varying amplitude and fixed exponent. Ra- 
dio quasars are splitted in bins of increasing radio luminosity 
[log(Pi.4GHz)>25.0,>25.5,>25.75]. Fits are calculated at two dif- 
ferent spatial scales, 0.1 < r p < 0.8 Mpc h~ 1 (solid) and 
1 < r p < 20 Mpc h~ 1 (dotted). A single fit over the range 
0.1 — 20 Mpc h _1 is indicated by the thick blue line, enclosed 
by the shaded error region. 



to probe. However, the clustering amplitude of RLQSOs 
increases as a function of radio luminosity, while that of 
RLAGNs decreases at the very highest radio luminosities. 
Extrapolation of our results suggests that both kinds of 
AGN might have similar clustering at radio luminosities 
in excess of 10 26 W Hz -1 . This would imply that the 
unified model for radio-loud quasars and radio galaxies can 
only be valid at these very high radio luminosities. This is 
consistent with the dependence of the relative numbers of 
the two AGN types as a function of radio power (Figure [Q. 



5 SUMMARY 

In this work, we have successfully applied cross-correlation 
techniques to characterize the environments of ~ 14, 000 
radio-loud AGN with Pi. 4 ghz > 10 24 W Hz^ 1 , selected from 
~ 1.2 million LRG at 0.4 < z < 0.8. We have also compared 
the clustering of RLAGN with that of radio-loud quasars 
over the same redshift interval. By using control samples of 
radio-quiet objects matched in redshift, stellar mass and op- 
tical luminosity (or radio luminosity, when appropriate) we 
have isolated the effect such parameters have in influencing 
the clustering signal. The main results of this paper can be 
summarized as follows: 

• Radio AGN at 0.4 < z < 0.8 are substantially more 
clustered than their parent luminous red galaxy popula- 
tion. Radio-loud AGN are also more strongly clustered than 
radio-quiet galaxies of the same stellar mass and redshift. 
The clustering differences are largest on scales less than 
1 Mpc ft -1 . 

• Radio-loud AGN hosted by more massive galaxies are 
more strongly clustered than those hosted by less massive 
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galaxies. However, the clustering difference between RLAGN 
and control samples of radio-quiet galaxies is most pro- 
nounced for RLAGN in low mass hosts. 

• We study the dependence of the clustering amplitude 
on the luminosity of the radio source. For r p > 1 Mpc 
there is a weak, but significant anti-correlation with radio 
power. For r p < 1 Mpc hT 1 the dependence of clustering am- 
plitude on luminosity is more complex: the cross-correlation 
amplitude increases with luminosity up to ~ 10 25 ' 3 W Hz -1 , 
and then decreases for the most luminous radio sources in 
our sample. 

• We have compared the environments of radio-loud AGN 
and radio-loud QSOs. RLAGN are clustered more strongly 
than RLQSOs on all scales, indicate that they populate dark 
matter halos of different mass. These results hold even when 
the RLAGN and RLQSO samples are matched in radio lu- 
minosity and black hole mass. 

• There are indications that the very most luminous 
RLAGN and RLQSOs in our sample (P> 10 26 W Hz" 1 ) do 
have similar clustering amplitudes. Only at these very high 
radio powers are the space-densities of radio-loud quasars 
and radio galaxies similar. This implies that unification 
of the two AGN populations can only be valid above P~ 
lO^WHz -1 . 

One major limitation of this study with regard to con- 
straining AGN unification scenarios, is that it is based purely 
on photometric data from the SDSS, so we are unable to split 
our RLAGN sample into high-excitation and low-excitation 
sources. It is quite possible that the presence or absence of 
emission lines will provide the best way to define a pop- 
ulation of radio galaxies that are clearly unified with the 
quasars. In this case, we would expect to find that the high- 
excitation radio galaxy population would cluster in a similar 
way to the quasars. 

In addition, we note that because the parent sample of 
our RLAGN catalogue consists of luminous red galaxies, it is 
also likely that we completely miss some number of RLAGN 
with bluer colors and stronger emission lines. The analysis of 
the R LAGN luminosity function presented in lDonoso et al.l 
(2009) indicates that the missing sources cannot constitute 
more than ~20% of the total RLAGN population, so will not 
dominate the clustering signal of the radio AGN population 
as a whole. Nevertheless the quasar analogues among the 
radio galaxy population may still be under-represented in 
our analysis. 

Fortunately, upcoming large spectroscopic surveys such 
as BOSS will target nearly complete samples of more than 
a million massive galaxies at 0.4 < z < 0.8 and will pro- 
vide optical spectra for tens of thousands of radio galaxies. 
We will then be able to quantify the fraction of RLAGN of 
given radio luminosity that have emission lines and how the 
clustering depends on emission line strength. 

The most definitive result to emerge from our analysis 
is clear proof that the environment of a galaxy on the scale 
of the dark matter halo in which it resides (i.e. on scales of 
~1 Mpc h -1 and below), does play a key role in determining 
not only the probability that a galaxy a is radio- loud AGN, 
but also the total lum inosity of the radio jet. Combining our 
results with those of iBest et al.l (|2005l ). we conclude that 
both black hole mass and environment must determine the 
radio-loud character of an active galaxy. 



Our previous work also demonstrated that strong evolu- 
tion of the radio AGN population only occurs above a char- 
acter istic radio luminosity of ~ 10 25 W Hz -1 (|Donoso et al.l 
2009). It is very intriguing that the results in this paper 
indicate that this luminosity marks the break point in clus- 
tering trends, and that the radio luminosity where denser 
environment ceases to have a boosting influence also is of 
order 10 25 W Hz" 1 . 

Finally, the strong evolution of the radio source popu- 
lation at radio luminosities above ~ 10 25 W Hz -1 combined 
with the strong clustering of this population, must imply 
that the heating rate of the gas in groups and clusters of 
galaxies is higher at redshifts ~ 0.5 than it is at the present 
day. We intend to quantify this in more detail in upcoming 
work. 
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